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Abstract. The aim of the paper is to derive the complex- valued least-squares 
estimator for bias-noise mean and variance. 



1. Kriging 

Let us consider a stationary random process e = {e^; N± 9 j D i = 1, . . . , n} with 
zero mean 

E{ (l } = E{ej} = E{e} = 
and the background trend J^ fe fjkP k = fjk/3 k (some known mean function m(j) 
with unknown regression parameters f3 k , where k = 1, 2, . . . ) then 

Vj = tj + m(j) = e.j + f jk /3 k e.g. j = n + 1 

and 

Vi = ei + m(i) = e,; + f lk /3 k i = l,... , n , 
where fjk (/') is a given vector and fik (F') is a given matrix. 
The unbiasedness constraint on the estimation statistics Vj = i^jV, 

E{V 3 } = E^Vl} 

produces the system of N(k) equations in the n unknowns 

fik = wj/ifc (/' = w'F) . 

For white noise 

" ej] 2 } = 0-2 + ^ 2 ^P«wj (£{[e - e] 2 } = <7 2 + ctVAw) , 

where pa (A) is the identity auto-correlation matrix, the minimization constraint 

dE{[e.j - tj] 2 } 2 2 fc 
— j = 2cr pjjWj + la fikfij = U , 

i 

where 

" ^] 2 } = °* + ^>;^ + 2a 2 (w*/ ifc - f jk ) p) , 

o 

let us add the n equations in the N(k) unknowns (ij to the system 

PiiUj) = -fikflj 

equivalent to 

^3 - -p JikHj 

substituting this term into the unbiased system 

fkj = fki^j 



2 



T. SUSLO 



WC get 

Pj = —{fkiP ll fik)~ 1 fkj 

and the kriging weights 

Uj = P U fikifkiP™ fik)^ 1 fkj = fjk{fkiP a fik)~ 1 fkiP lt ■ 

Now, we can write the kriging estimator 

v 3 = u) Vi = f 3k p k (v = w'v = /'/?) , 
where the least-squares estimator 

P k = {f kl p ll f l k)- 1 fk lP ll v t (p = {F'k^F^F'k^v) 
is the best linear unbiased estimator for n — > oo 

Urn E{[Vj - f jk (3 k } 2 } = lim E{[v d - f jk (3 k } 2 } = , 

where 

fjkP k ] 2 } = O 2 0j)piiU) = -a 2 u)Ukp) = -° 2 !-jkp) = & 2 fjk (fkiP™ fik fkj ■ 

2. Complex-valued bias-noise mean and variance 
Since for constant bias-noise mean (k = 1) 

Vj = ej + f kj (3 k = € d + P 1 

and 

Vi = e, + f kl f3 k = a + (3 1 i = l,...,n 
the precession of the estimation statistics can not be compared to zero value 

Envj-a 1 } 2 } 

let us introduce the bias-noise mean with non-zero slope (k = 1,2) 

V j =e j + f jk f3 k = e j +[3 1 +jf3 2 

and 

Vi = e, + f ik /3 k = a + I3 1 + i(3 2 i = l,...,n 
to find the best linear unbiased estimator for any n > 2 

E{[Vj f lk f3 k } 2 } = Eilvj f jk (3 k ] 2 } = 

we have to fulfill 



a 

n 



fjkifkiP™ fik) 1 fkj — [ 1 j 



ni ni z 



f - 2m n j + m s 



naf, 



at 



in 



j =m n ± Icr„ 



v, = w)vi = f 3k p k - /3 1 + fp 2 = b + ja, 

where: I = y^T, rn n = i = ± £\ i, m sn = i 2 = \ J^i °n = \J i 2 - and we 
get simple mean 



in 



= b + m n a ± Ia n a = Vi ± Ia n 1 (ivi - m n Vi) 
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and variance 

-2 i 2 - 2 

a = uj v i - m , 

where 

Lu l vf = vf ± la' 1 (ivf - m n vf) , 
charged by the imaginary error. 

Example. Let us consider the best linear unbiased estimator of mean 

W=»(m) = 6.1 

for a test sample yi{xi) 



Xi 


1.7 


2.1 


3.9 


7.2 


8.6 


8.5 


7.3 


5.1 


2.8 


1.8 


1.6 


Vi 


3.2 


3.9 


4.9 


5.3 


5.5 


6.2 


6.5 


6.9 


7.5 


8.3 


9.4 



with the real-valued standard error 



±\ M =±0.5 
V n 

and the imaginary one 

3?(m) = ±a^/x? - xf = TO. 2 . 
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